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摘 要 : 针对 生 鲜 闭环 供应 链 网 络 设计 问题 ， 建 立 了 一 种 基于 生 鲜 闭环 供应 链 网 络 的 便 棒 优化 模型 ， 以 解决 供应 
链 网 络 中 的 不 确定 性 问题 。 首 先 ， 针 对 涵盖 5 个 节点 的 生 鲜 供应 链 网 络 结构 ， 建 立 了 多 周期 ， 多 产品 ， 以 最 小 化 成 
本 ， 最 小 环境 影响 为 目标 的 混合 整数 规划 模型 ， 采 用 模糊 折 中 规划 与 区 间 数 据 鲁 棒 优化 方法 进行 处 理 ; 其 次 ， 在 
原 有 蜜 儿 算 法 的 基础 上 ， 引 入 差分 进化 原则 ， 增 强 算法 的 全 局 搜索 能 力 与 收敛 速度 ; 最 后 ， 通 过 Matlab 数值 分 析 
与 仿真 实例 表明 ， 所 提 章 棒 优 化 模型 与 密 获 算法 在 求解 生 鲜 闭环 供应 链 网 络 设计 问题 中 具有 明显 优势 。 
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Robust optimization design of fresh closed-loop supply chain network based on 
improved honey badger algorithm 


Dong Haia, Lin Guodong® 
(a. School of Applied Technology, b. School of Mechanical Engineering, Shenyang University, Shenyang 110044, China) 


Abstract: Aiming at the design of the fresh food closed-loop supply chain network., this paper introduced a robust optimization 
model based on the fresh food closed-loop supply chain network, which was established to solve the uncertainty problem in 
the supply chain network design. Firstly, for the network structure of fresh food supply chain covering 5 nodes, with the goal 
of minimizing cost and minimizing environmental impact, a mixed integer programming model is established, and the 
compromise planning and interval data robust optimization method are used for processing. Secondly, based on the original 
honey badger algorithm, the principle of differential evolution is introduced to enhance the global search ability and 
convergence speed of the algorithm. Finally, it is shown that there are obvious advantages for the robust optimization model 
proposed in this paper and the differential honey badger algorithm can solve the closed loop of fresh food in supply chain 
network design problems efficiently through Matlab numerical analysis and simulation examples. 
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0 引言 的 生 鲜 电 商 供应 链 决 策 模 型 ， 并 设计 “收益 共享 -双向 成 本 

分 担 ” 契 约 ， 实 现 了 生 鲜 电 商 供应 链 的 完美 协调 与 帕 紧 托 改 

由 于 自然 资源 短缺 、 工 业 产 品 残留 对 人 类 生活 和 环境 危 。 进 。 王 椅 等 门 构建 了 受 生 鲜 农产品 新 鲜 度 和 价格 影响 的 消费 

二 害 等 问题 的 出 现 ， 逆 向 物流 和 闭环 供应 链 的 设计 已 引起 研究 ”者 时 变 效 用 函数 ， 并 通过 建立 由 零售 商 和 供应 商 组 成 的 两 级 

三 ”人员 和 决策 者 的 高 度 重视 。 生 鲜 供应 链 因 其 易 损耗 ， 难 调节 生 鲜 农产品 供应 链 利润 模型 ， 实 现 供应 商 的 最 优 保鲜 努力 和 
@ 等 问题 逐步 成 为 供应 链 领 域 的 研究 热点 。 零售 商 的 最 优 定价 。 

国内 外 学 者 对 于 生 鲜 供应 链 的 研究 方向 大 多 集中 在 定价 在 闭环 供应 链 网 络 研 究 方面 ， 孙 嘉 轶 等 四 基于 消费 者 偏 

与 库存 策略 、 保 鲜 努 力 和 协调 优化 等 方面 。 冯 颖 等 由 在 考虑 ”好 与 公平 关切 ， 建 立 了 单 双 渠 道 的 闭环 供应 链 网 络 决策 模型 ; 

生 鲜 农产品 随机 产 出 和 价值 损耗 等 特性 前 提 下 ， 分 别 构建 了 张 诸 等 四 针对 不 确定 条 件 下 的 可 持续 闭环 供应 链 网 络 设计 的 

离 岸 价格 和 到 岸 价 格 两 种 模式 下 的 分 散 决 策 博弈 模型 。 马 雪 ”问题 ， 建 立 了 以 最 小 成 本 和 环境 影响 、 最 大 社会 影响 为 目标 


丽 等 包 研 究 三 级 供应 链 冷 链 系 统 下 由 第 三 方 物流 承担 生 鲜 品 。 且 带 有 模糊 参数 的 多 目标 闭环 供应 链 网 络 规划 模型 ， 董 海 等 
的 保鲜 与 低 碳 责任 时 ， 冷 链 系 统 的 保鲜 、 碳 减 排 及 定价 等 19 研究 了 电网 中 断 条 件 下 的 生 鲜 闭环 供应 链 网 络 多 目标 模 
策 问题 ， 提 出 了 批发 价格 和 两 部 收费 相 结 合 的 契约 激励 机 制 。 糊 优 化 设计 模型 ,解决 了 生 鲜 供应 链 网 络 设计 中 的 不 确定 性 
Mohammadi 等 上 针对 生 鲜 产品 易 腐败 ， 浪 费 等 问题 ， 开 发 问题， 朱 晨 等 5 研究 了 由 零售 商 主 导 的 双 渠 道 闭 环 供应 链 
了 一 种 基于 保存 技术 投资 的 新 型 协调 机 制 ， 提 高 了 整体 供应 ”协调 问题 ， 采 用 利润 共享 一 费用 分 担 契 约 实现 供应 链 的 协调 ; 
链 的 利润 水 平 。 闻 弄 等 和 由 针对 自然 灾害 对 产 出 的 影响 因素 ， 王道 平等 (4 研究 了 信息 非 对 称 下 考虑 制造 商 回收 行为 的 闭 
建立 基于 政府 、 零 售 商 和 农户 的 三 阶段 Stackelberg 博弈 模型 ，“ 环 供应 链 协调 问题 ， 运 用 伊 蕨 过 程 阐 释 回收 率 的 随机 演化 过 
对 比分 析 在 政府 不 同 补贴 政策 以 及 零售 商 的 合作 偏好 行为 下 程 ， 构 建 分 散 决策 模型 ， 获 得 制造 商 与 零售 商 的 最 优 均 衡 解 ; 
政府 的 最 优 补贴 率 和 零售 商 的 最 优 收 购 价 格 。 邱 慧聪 等 加 针 ”Yavari 等 03] 研 究 不 确定 条 件 下 易 腐 产 品 的 绿色 闭环 供应 链 网 
对 碳 税 政策 的 生 鲜 绿色 供应 链 网 络 优化 问题 ， 设 计 了 成 本 与 ” 络 的 设计 问题 ， 设 计 了 创新 型 混合 整数 规划 和 鲁 棒 优化 模型 ， 
碳 排放 量 平衡 优化 的 生 鲜 品 绿色 供应 链 网 络 结构 ， 并 建立 多 ”Dey 等 4 分析 了 在 闭环 供应 链 网 络 设计 中 制造 商 作 为 斯 塔 克 
目标 混合 整数 模型 。 刘 墨 林 等 四 通过 构建 集中 式 与 分 散 式 下 尔 伯 格 领导 者 时 ， 两 个 竞争 零售 商 的 不 同 博弈 策略 。 
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规模 优化 问题 Pa， 
大 规模 问题 时 


or 


定稿 


闭环 供应 链 网 络 


。 近 年 来 ， 该 领域 的 最 新 下 
生 优 化 问题 PD 7。 


搜索 行为 ， 因 其 实验 结 


良好， 结构 简章 


算法 在 处 型 


泛 且 良 好 的 应 用 前 景 。 
局 搜索 能 力 ， 收 敛 速度 等 方 


下 还 有 一 定 的 改进 
高 该 算法 的 收敛 速度 和 


链 网 络 多 


综 上 所 述 ， 生 鲜 


久 献 相对 较 少 ， 且 对 生 


的 不 确定 诉 
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题 是 一 个 典型 的 NP-hard 问 
多 采用 智能 算法 来 求解 此 类 
法 (HBA) 是 由 Hashima 等 09] 于 
年 提出 。 该 算法 模拟 了 蜜 获 的 挖掘 和 寻找 蜂蜜 的 动态 
名， 可 以 有 效 求解 大 


然而 原始 


空间 ， 可 将 其 与 其 他 算法 或 理念 相 结合 ， 
1 寻 优 能 力 。 
闭环 供应 链 网 络 优化 设计 的 相关 研究 和 
FE 参 数 的 处 理 方式 略 


单一 。 因 此 ， 本 文 将 建立 一 个 考虑 环境 影响 的 生 鲜 闭环 供 


换 为 单 目标 模型 ， 再 通过 区 间 重 棒 亿 


数 进行 处 理 ; 在 原 
行 改进 以 求解 所 提 模 型 ， 最 后 通过 实例 验 订 
的 可 行 性 与 算法 的 优越 性 。 


对 其 进 


计 问 题 。 其 化 


模型 建立 


问题 描述 


标 优 化 模型 ， 并 采 


基础 


本 文 将 研究 由 供应 商 、4 
心 组 成 的 多 产品 、 多 时 期 、 多 梯队 的 生 鲜 闭 
应 链 网 络 运 行 如 下 : 首先 通过 


络 与 运输 方式 ， 将 4 
产品 通过 仓库 从 制造 商 运送 到 零 外 


E 产 所 需 的 材料 从 供应 


E 产 商 、 仓 库 、 零 售 商 和 收集 中 
] 环 供应 链 网 络 设 
见 有 的 供应 链 网 
商 运送 到 生产 商 ; 


模糊 规划 方法 将 多 目标 转 
化 对 模型 中 的 不 确定 性 
上 上 ， 引 入 差分 进化 原则 
E 所 提 模 型 
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参数 都 是 预 


保存 在 收集 中 心 。 在 收集 中 心 ， 


和 商 。 不 满意 的 产品 被 退 世 
将 可 用 的 退回 物品 作为 半 


成 品 发 送 给 制造 商 。 其 


商 的 需求 与 退货 率 处 于 一 个 


不 确定 性 的 状态 。 有 具体 供应 链 网 络 结构 如 图 1 所 示 。 
画 一 > 
百 一 > 虚线 为 逆向 物流 
此 = et 实 线 为 正 向 物流 
供应 商 8 制造 商 M 仓库 W 
a a 
a gE 
一 一 多 © 一 > 
m 一 
di .mm 
上 零售 商 R 客户 群体 


图 1 生 鲜 闭环 供应 链 网 络 结构 示意 图 


Fig.1 Schematic diagram of fresh closed-loop 


supply chain network structure 


1.2 模型 假设 及 符号 定义 


本 文 假设 如 下 : a) 所 有 设施 的 潜在 位 置 、 容 量 和 成 本 


E 确 定 的 ，b) 每 个 零 


心 ， 每 个 零售 商 


有 需求 ;，c) 零售 商 退 


售 商 仅 从 一 个 仓库 接收 


有 产品 只 能 送 到 


每 个 产品 的 百分比 作为 退回 


只 将 


单位 运输 成 本 是 固 
从 仓库 到 零售 商 以 及 从 零 


应 商 到 制造 商 、 从 制造 商 到 


定 的 ， 


P 心 之 间 的 每 个 时 期 都 


品 发 送 回收 集中 心 ; 
仓库 以 及 从 收集 中 心 到 
其 运输 成 本 与 运输 量 了 


个 收集 中 


的 产品 发 送 到 一 个 收集 中 心 ; dd) 


e) 从 供 
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参数 定义 如 下 : 五 表示 产品 i 的 退货 紊 ; L44 表 示 收 集中 
心 a 的 容量 ; LMm 表示 制造 商 m 的 生产 能 力 ; LS; 表示 供应 
商 s 的 供 货 能 力 ; LW 表示 仓库 w 的 容量 ，pi 表示 产品 i 的 
生命 周期 ，FC44a、FCMm 以 及 FCW 分 别 表 示 建 立 收集 中 心 
a、 制 造 商 m、 仓 库 w 的 固定 成 本 ; COMimi 表 示 制 造 商 m 在 
t 时 间 段 内 生产 产品 i 的 单位 成 本 ; DR 表示 零售 商 r 在 1 时 
间 段 对 产品 i 的 需求 ， E44、EMm、EW, 分 别 表示 建立 收 负 
中 心 a， 制 造 商 m 和 仓库 w 的 环境 影响 ;El4ia、EIMim、 
EIRi、EIWiw 分 别 表示 产品 i 储存 在 收集 中 心 a， 制 造 商 m， 
零售 商 " 以 及 仓库 w 四 处 的 环境 影响 ，ET4jam 表示 将 产品 j 
从 收集 中 心 a 运往 制造 商 m 的 环境 影响 ，ETMimw 表示 将 产 
品 i 从 制造 商 m 运输 至 仓库 w 的 环境 影响 ，ETSsm 表示 供应 
商 s 向 制造 商 m 运输 的 环境 影响 ，ETRAs 表 示 预 分 配 卡车 从 
零售 商 " 到 收集 中 心 a 行程 的 环境 影响 ，ETW 表示 预 分 配 
的 卡车 从 仓库 w 到 零售 商 r 行 程 的 环境 影响 ，HC4ia、HCMim、 
HCRi、HCWiw 分 别 表示 产品 i 在 收集 中 心 a、 制 造 商 m、 零 
售 商 r 以 及 仓库 w 四 处 的 单位 库存 成 本 。 

Oi 表示 产品 i 的 使 用 系数 ， 若 产品 j 能 作为 产品 i 的 原 
材料 ， 则 0; 为 1， 否则 为 0;，Goi 表示 生产 其 他 产品 所 需 的 退 
货 产 品 i 的 百分比 PCRAia、PCSMsmt 表示 在 t 时间 内 ， 为 
中 心 a 从 零售 商 r 处 购买 产品 i 的 单位 采购 成 本 以 及 为 
制造 商 m 从 供应 商 s 处 采购 原材料 的 单位 采购 成 本 ; 
TCAMiam、TCMWimwm 分 别 表 示 在 t 时 间 内 ， 产 品 i 从 收集 中 
心 a 运送 至 制造 商 m 以 及 从 制造 商 m 运 至 仓库 w 的 单位 运 
输 成 本 ; TCSsmi 表示 在 t 时 间 内 ， 从 供应 商 s 运送 至 制造 商 
m 的 原材料 的 单位 运输 成 本 ; TCRAwa、TCWRw 分 别 表示 在 
1 时间 内 ， 从 零售 商 r 到 收集 中 心 a 以 及 从 仓库 w 到 零售 商 
r 的 固定 运输 成 本 ; 若 产 品 能 从 零售 商 x 运送 至 仓库 w， 则 
Um 值 为 1， 反 之 为 0; 若 产品 能 从 零售 商 r 运送 至 收集 中 心 
a， 则 Vs 值 为 1， 否 则 为 0。 
决策 变量 符号 定义 如 下 : Nimi 表示 在 t 时 间 内 ， 制 造 商 
m 生产 的 产品 i 的 数量 ， IL4iar、ILMimr、ILRim、ILWm 分 别 
表示 在 1 时间 段 内 ， 产 品 i 在 收集 中 心 a、 制 造 商 m、 零 售 商 
r 以 及 仓库 w 中 的 库存 水 平 ，Nsmi 表示 在 t 时 间 内 ， 从 供应 
商 s 运 至 制造 商 m 的 原材料 数量 ; NWRiwn、NRAirat、NAMiamt、 
NMWimwi 分 别 表示 在 t 时 间 内 ， 产 品 i 从 仓库 w 发 货 至 零售 
商 x、 从 零售 商 r 运送 到 收集 中 心 a、 收 集中 心 a 运送 到 制造 
商 m 以 及 从 制造 商 m 运 至 仓库 w 的 产品 数量 ;O04,OM,OW 
以 及 OR 分 别 表示 收集 中 心 ， 生 产 商 ， 仓 库 以 及 零售 商 的 状 
态 ， 若 上 述 四 个 供应 链 节 点 处 于 开放 状态 则 值 为 1， 否 则 为 
0。 目 标 函 数 Z1 为 最 小 化 网 络 成 本 ，Z 有 为 最 小 化 环境 影响 (以 
CO 排放 量 为 衡量 标准 )。 

ZI= 各 设施 的 总 库存 成 本 + 建立 制造 商 、 仓 库 和 收集 中 
心 的 总 固定 成 本 + 设施 之 间 的 总 运输 成 本 + 生产 加 工 各 类 产 


7 


污 
浪 


制造 商 的 


商 到 收集 中 心 的 运输 成 本 是 静态 


的 ， 与 运输 量 无 关 ， 预 先 分 配 的 卡车 在 仓库 、 零 售 商 和 收集 
行程 ，g) 零售 商 处 于 固定 位 置 日 


设施 的 能 力 ( 供 应 商 、4 
限 的 。 


E 产 商 、 仓 库 、 零 售 商 和 收 


集 中 心 ) 是 


本 文 以 乳 制品 供应 链 网 络 为 例 ， 定 义 符号 ， 相 关 参 数 、 


和 决策 变量 如 下 所 示 。 
的 潜在 位 置 的 集合 ， 


付 写 契 


义 如 下 。4 代表 收集 中 心 
aE4; I 代表 产品 ，iE7; 世 和 产 分 别 代 
表单 源 产品 和 再 利用 产品 ，is EIT，ic ET; M 代表 制造 商 的 潜 


集合 ，m EM; R 表示 零售 商 的 固 


3 表示 供应 商 的 回 定 地 点 得 


表示 仓库 的 潜在 位 置 


定 地 点 集合 ， rEeR; 
晨 合 ，s ES; 1 表示 时 间 段 ,tfE7; WW 
， WEW, 


品 的 总 生产 成 本 。 
DZ= 各 设施 保持 库存 时 的 CO> 排 放量 + 建 立 制造 商 、 仓 库 
和 收集 中 心 的 CO2 排 放量 + 运输 的 CO 排放 量 。 
HinD= (YY HCM;, x 1M + >》FCIxEWv+> > YHCR, x 


ILRa + >,) YHCA, x LA + (2, FCM,, x OM + ,FCW x 
"J. 三 m w 


YD Drm, Ox NAM jn) + (1) 
Fa Wr 
FE PR x NRAsa + 5 YD PCSM, x Ns)+ 
区 P a 1 党 m 
(六 >》》》03xCOMiw x(》 QO; x NAM jon) + 
和 a m 1 
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> > 2 EA XILA + > 2 2 EIR, xILR,)+ 
(LEM, xOM + >,EW, xOW+2,EA, xOA)+ 
CC > > ETS,, x Ns + 2 > 2 ETM, XNM Wi + O) 


DETW,, xU,, + > > ETRA, x Va t+ 


加 


式 G) 和 (人 ) 分 别 表示 零售 商 和 仓库 中 有 充足 的 输入 。 
二 DR < ENWR, ,Vr G) 
TEN > ENR Win @) 


式 (5) 表 示 从 每 个 零售 商 返 回 到 收集 中 心 的 产品 数量 都 
是 其 需求 的 一 部 分 。 


DD NRA < > DR xk Vir (5) 
式 (6) 表 示 每 个 收集 中 心 的 输入 流 大 于 输出 流 。 
> > axMNR4w > NAMi, via (6) 
式 (7) 表 示 每 个 制造 商 的 输入 和 输出 之 间 的 平衡 。 


2 Nm > 2 2 20x NMWom -2 2 20 * Ox NAM mm (7) 
t 下 并 


学 


式 (8) 讨 论 了 每 个 制造 商 的 产量 至 少 应 该 与 从 制造 商 送 
到 仓库 的 产品 一 样 多 。 


之 Ni > > 之 NMWw sism (8) 


式 (9) 至 (12) 分 别 表示 制造 商 、 供 应 商 、 仓 库 和 收集 中 心 
的 容量 限制 。 


> Ni < IM, x OM ,Vmt (9) 

SN <LS,,Vs,t (10) 

We + 2 2 NMWom <LW, xOW, Vw,t (11) 
2 Aa + 2 2 NRAw <LA, x OA, Va,t (12) 


式 (13)~(16) 表 示 制 造 商 、 仓 库 、 
间 段 内 每 种 产品 的 库存 水 平 。 


收 货 中心 、 零 售 商 各 时 


TM = Nog t ILM sot) = ZS NMWom Vism,t (13) 
ILW,, = 之 NM 二 下 了 一 之 NIWR Vi (14) 
ILA, = NRA +ILA -2 NAMi , Vi,a,t (15) 
ILAia = > NWR, + LR 一 DR .Vis r,t (16) 
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式 (17) 确 保 产 品 只 能 由 现 有 的 制造 商 生 产 。 


OM £2 Nm < MxOM, Vm (17) 
式 (18)(19) 确 保 产 品 从 生产 商 发 送 到 仓库 或 收集 中 心 。 
OW 222 NMWo SM xOM,vw (18) 
04 < YY NRA < M04, va (19) 
式 (20)(21) 确 保 任何 非 链接 设施 之 间 没 有 转运 物流 。 
U,, < 2 NWR <M xU, ,Vw,r (20) 
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Vs<s> >》NR4Aw < M xV,, ,Vr,a CD 


式 (22) 表 示 每 个 零售 商 只 从 一 个 仓库 接收 产品 ， 式 (23) 
表示 每 个 零售 商 只 将 返回 的 产品 发 送 到 一 个 收集 中 心 。 


YU <1, Wr (22) 


之 多 Ls VF (23) 


式 (24) 表 示 收 集中 心 发 送 的 所 
间 被 用 于 生产 再 制造 产品 。 
NAM > >》 YQ; x NAMi » Vim,t 024) 


退货 产品 必须 在 同一 时 


式 (25) 表 示 产 品 i 在 每 个 仓库 的 库存 水 平 在 下 一 个 Px( 产 
品 i 的 生命 周期 ) 连 续 期 间 内 始终 小 于 产品 i 的 总 出 库 量 。 


t+(pi-l) 
ILWi < > NWRiw , Visw,t (25) 
KE 二 


式 (26) 确 保 零售 商 在 每 个 时 期 的 每 个 产品 的 库存 水 平 小 
于 下 一 个 Pi 连续 时 期 的 总 需求 。 


t+(pi-l) 
ILR, < > DR Visr,t (26) 
$ 


2 ”模型 处 理 


2.1 多 目标 处 理 
由 于 本 文 所 建 的 两 个 目标 函数 1 与 甩 在 某 种 程度 上 存 
在 冲突 ， 即 2Z1 目标 函数 表现 较 优 时 可 能 会 导致 男 外 的 子 目 
标 函 数 之 表现 不 佳 。 为 了 实现 帕 累 托 最 优 ， 本 文 将 采用 模 
糊 折 中 规划 方法 进行 多 目标 处 理 。 根 据 文献 [20] 中 描述 ， 模 
糊 折 中 规划 法 的 具体 数学 表达 式 如 下 : 


w|i 07) 


其 中 p 属于 [1,%] 中 的 正 整 数 ， 良 表示 各 个 子 目 标 函 数 记 的 对 
应 的 理想 解 ， 即 =min(f); 大 表示 赋予 各 个 子 目 标的 权重 ， 
且 大 属于 (0，o)。 从 而 本 文 的 多 目标 模型 将 等 价 转换 如 下 : 
mn] (2 
其 中 wi 为 各 个 子 目 标 函数 的 相对 权重 ， 且 好 为 友 与 六 的 乘 
积 。 通 过 改变 参数 p 的 值 可 以 得 到 不 同 的 有 效 解 ， 最 常见 的 
值 是 p=1、2 和 = 。 
2.2 鲁 棒 处 理 

为 了 将 确定 性 模型 扩展 到 不 确定 环境 中 ， 本 文 假设 退货 
率 、 需 求 为 不 确定 性 参数 ， 鲁 棒 方法 主要 有 两 阶段 : 首先 是 
构建 不 确定 参数 区 间 对 不 确定 参数 进行 处 理 描述 ， 其 次 是 采 
j 鲁 棒 对 等 式 进行 转换 ， 将 鲁 棒 模 型 转换 为 等 价 的 混合 整数 
规划 模型 进行 求解 。 
本 文采 用 设置 不 确定 集 区 间 考 虑 不 确定 性 参数 的 偏差 范 
所 RI， 提出 闭环 供应 链 模 型 的 鲁 棒 对 应 形式 。 此 外 ， 本 文 
定义 了 参数 d=DRit，k=Ki。 不 确定 参数 DRi，Ki 具有 相同 
独立 分 布 ， 并 且 在 范围 [4-44+4]，[*-ht+4] 中 分 布 对 称 。 此 
外 ， 不 确定 性 预算 了 取 值 [0,1]。 
于 假设 需求 与 退货 率 是 不 确定 的 ， 并 且 考 虑 鲁 棒 优化 
模型 中 最 坏 情况 参数 对 模型 和 算法 的 影响 ， 即 本 文 只 考虑 需求 
和 退货 率 的 不 确定 参数 的 正 偏差 。 不 确定 性 参数 主要 集中 于 约 
束 式 (3)(5)(26)。 以 约束 式 (3) 为 例 ， 该 约束 左 侧 含有 不 确定 参 
数 ， 将 其 不 确定 参数 值 记 为 4;, jEJ， 而 /表示 这 个 约束 中 不 
确定 性 参数 的 数量 索引 集合 ， 根 据 文献 [2H] 定 义 4=4+45 ， 
其 中 尹 被 称 为 标 称 值 ，4 表示 4 的 最 大 偏离 值 ， 5 EU，U 表 
示 一 个 有 界 的 闭 集 ， 而 性 则 在 不 确定 集 U 中 选择 以 控制 少 的 
上 下 限度 。 根 据 上 述 定 义 内 容 ， 约 束 式 (3) 可 表示 为 
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dx 这 pi DINWR, ,VE 


其 中 : jy 表示 需求 的 保守 性 参数 ,ju Ee [0,lilx|r|]。 

综 上 所 述 ， 引 入 4 表示 供应 链 中 产品 需求 的 不 确 
即 a=4dj+dt)， 和 EU，jEJ。g 为 控制 @ 的 独立 随机 
根据 上 述 文献 定义 及 规则 ， 引 入 辅助 变量 六 。 此 外 ， 
性 预算 用 取 值 [0,1]。 
可 得 到 相关 约束 与 鲁 棒 等 价 模型 


ES a 1+(pi-1) 
e=DRmxKi; B= >》 DR，。 
1 


Z* 如 下 。 


约束 中 


等 : 基于 改进 HBA 算法 的 生 鲜 闭环 供应 链 网 络 鲁 棒 优化 设计 


(29) 


定 中 ? 
变量 。 
不 确定 


除 此 之 外 ，ME [0,lilx|r|x| 相 。 
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式 ” 和 “蜂蜜 模式 ”该 算法 使 用 改变 搜索 方向 的 标志 了， 
确保 能 有 更 多 的 机 会 对 搜索 空间 进行 严格 扫描 ， 以 跳出 局 章 
最 优 的 困 局 。 

1) 挖 掘 模式 

在 挖掘 模式 中 ， 窒 多 的 运动 路 线 类 似 于 心 形 ， 通 过 式 


G7) 可 以 模拟 运动 路 线 轨迹 。 
prey + PX PXTX Xe 十 
Fxnxaxd;xlcos(2xn)x[1—cos(2xr)] 
Oy xnew 代表 蜜 获 的 新 位 置 ，xprey 代表 猎物 的 位 置 全 局 最 
佳 位 置 ，f21( 默 认 值 为 6) 描述 的 是 蜜 获 获得 食物 的 能 力 ， 


Nnew = 


(G37) 


Min{Z1, 22} di 是 猎物 和 第 i 只 蜜 获 之 间 的 距离 ，r3, m4 和 rs 是 0 到 1 之 间 
PE G0) 的 三 个 不 同 的 随机 数 ，Pi 作为 改变 搜索 方向 的 标志 ， 由 式 
所 | 气 (38) 确 定 : 
Ee <o1+exi Es) Vi G1) hal a (38) 
加 re 属于 0 到 1 的 随机 数 ， 蜜 牙 在 挖掘 模式 中 主要 依赖 于 
DD ss 6 对 猎物 的 嗅觉 强度 矿 猎物 的 距离 d 和 时 间 搜索 影响 因子 a 
0<Y, <d;0<¢,<J (33) 进行 搜索 。 
其 他 约束 条 件 不 变 。 2) 蜂 蜜 模式 
条 个 PA 有 导 密 狂 、 梨 的 情 ; 本 | 
3 ”求解 方法 | 蜂蜜 模式 中 ， 蜜 效 跟 随 导 蜜 鸟 到 达 蜂 梨 的 情况 可 以 模拟 
为 式 (39)。 
3.1 编码 处 理 Now = Xprey + XH XAXd; (39) 
考虑 到 供应 链 网 络 设 计 问 题 中 的 多 层级 、 多 周期 的 复杂 其 中 ， 六 为 0 到 1 的 随机 数 ，x 和 局 由 式 (36)(38) 确 定 。 
性 ， 本 文选 用 矩阵 实数 编码 方法 : 假设 种 群 规模 为 P'"， 第 0 式 (39) 可 以 看 出 ， 蜜 牙根 据 距离 信息 di:， 在 目前 发 现 的 猎物 
代 种 群 克 o={Ki,K2,.….,KPp"}， 其 中 Kj 代表 第 0 代 的 第 j 个 个 位 置 xprey 附近 进行 搜索 。 在 这 一 阶段 ， 搜 索 行为 会 受到 随时 
体 ，jE[l P]，QE[l,tmax]; 定义 区 =(Cap)nxn，Cab 表 示 和 矩阵 ” 间 变 化 因子 (@) 的 影响 。 
元 素 ,，aE[1,m]; bE[1,n]。 除 此 之 外 ， 和 矩阵 万 的 大 小 跟 供 3.3 ”差分 富 获 算法 
应 链 网 络 中 的 层级 、 节 点 设置 数 有 关 ， 日 矩阵 元 素 Cu 主 为 改进 原始 蜜 获 算法 处 理 大 规模 问题 时 易 陷入 局 部 最 优 ， 
要 反映 与 上 层 供应 者 的 物流 或 者 供 货 关 系 。 例 如 到; 为 一 个 3 收敛 速度 慢 等 缺陷 ， 本 文 将 差分 进化 算法 (Differential 


X4 阶 的 和 矩阵， 则 表示 该 供应 链 网 络 模型 分 为 3 
多 设 
别 来 自 上 层 供应 者 第 1 个 节点 53 的 单位 供 货 量 与 第 3 
48 的 单位 供 货 量 ，0 则 表示 该 两 节点 间 没 有 供 货 关系 。 
3.2 ”原始 宣 狼 算法 


民 ， 节 点 数 最 
为 4，C21=(53,0,48,0) 则 表示 第 2 层 的 第 1 个 节点 接收 分 


个 节点 


Evolution A DEA) 中 的 交叉 变异 原则 引入 到 原始 蜜 
获 算法 中 ， 形 成 差分 蜜 获 算法 (Differential Honey Badger 
Algorithm，DHBA)。 在 选择 新 位 置 与 个 体 时 ， 通 过 执行 差 
分 进化 交叉 变异 策略 ， 使 得 DHBA 具有 全 局 搜索 范围 更 广 、 
种 群 更 具 多 样 性 。 具 体 DE/rand/1/bin 交叉 变异 策略 如 下 : 


3.2.1 种 群 的 初始 化 1) 变 异 策 略 
HBA 的 基本 思想 是 通过 模拟 蜜 获 的 疯 食 行为 进行 优化 。 Vpj = Xa + Px rand Xs —X9,)) (40) 
式 (34) 初 始 化 蜜 牙 的 数量 (种 群 大 小 N) 及 其 各 自 的 位 置 。 其 中 ， 训 ey 代表 当 六 前 搜索 区 域内 的 最 优 位 置 ， 下 为 缩放 系数 
X= lbi+ rix (ubi—lb) (34) ”上 且 为 (0，]) 间 的 常数 ， 有 sj、ZY9,; 为 种 群 中 的 随机 个 体 ，rs、 
其 中 ; zi 为 [0,1] 的 随机 数 。xi 代表 种 群 N 中 第 i 个 蜜 获 的 位 。 ro 属于 种 群 N 中 的 互 不 相同 的 整数 。 
; 而 1b: 和 upbi 分 别 是 搜索 域 的 下 限 和 上 限 。 2) 交 叉 操 作 
3.2.2 捕食 强度 小 py» if rand(0,1) < CR G41) 
蜜 获 的 捕食 强度 与 猎物 的 集中 强度 以 及 其 与 第 i 个 蜜 牙 X aj» otherwise 
的 距离 有 关 。 五 是 猎物 的 气味 强度 ， 如 果 气 味 高 ， 运 动 就 快 ， 其 中 rand(0,1) 是 (0,1) 服 从 上 均匀 分 布 的 随机 数 ，CR 为 交叉 
反之 亦 然 。 具 体 定义 如 下 : 概率 且 属 于 [0,1]。DHBA 具体 算法 流程 如 图 2 所 示 。 
= G5) 4 算 例 验证 
其 中 : 5S 代 表 猫 物 的 集中 强度 与 气味 强度 ，di 表 示 猎 物 和 第 i 只 本 文 以 沈阳 某 奶 制品 制造 厂商 为 例 ， 该 公司 供应 链 网 络 
猿 之 间 的 距离 ，S=( 人 -5 4=xwo-s。 7 为 [0.1] 的 随机 数 。 范围 主要 履 盖 市 内 五 个 行政 区 ， 由 5 个 节点 组 成 ， 其 中 包括 
3.2.3 更 新 密度 因子 3 个 供应 商 、5 个 制造 中 心 、4 个 仓库 、17 个 零售 商 以 及 3 
密度 因子 (oa) 控制 时 间 变 化 的 随机 性 ， 以 确保 算法 两 阶 ”个 收集 中 心 。 该 公司 主要 生产 两 类 奶 制品 ， 回 收 中 心 只 回收 
段 的 平稳 过 渡 。 根 据 式 (36) 更 新 随 迭 代 次 数 减 少 的 递减 因子 其 中 一 类 产品 ， 且 只 考虑 固定 单一 的 运输 方式 。 除 此 之 外 ， 
a， 减 少时 间 的 随机 性 。 根据 企业 实际 情况 ， 各 个 参数 值 都 服从 均匀 随机 分 布 ， 需 求 
scxom (> G0) 不 确定 性 DRimE [400,2000]， 退 货 率 不 确定 性 Kie[0,1]， 供 
fw 应 链 网 络 具 体 参数 值 设置 如 下 表 1 所 示 。 


其 中 : 
认 值 设 为 2)。 

3.2.4 更 新 粒子 位 置 
正如 上 文 所 述 ，HBA 粒子 位 


更 新 过 程 采 用 “ 挖 


fmax 表示 最 大 欠 代 次 数 ; C 表示 一 个 大 于 1 的 常数 ( 默 


掘 模 


4.1 和 鲁 棒 模型 验证 

为 验证 本 文 所 设计 的 鲁 棒 优 化 模型 的 稳定 性 与 可 行 性 ， 
针对 需求 不 确定 的 情况 生成 了 表 2 中 的 5 个 测试 问题 ， 相 应 
参数 如 表 1 所 示 。 
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供应 链 相关 参数 设置 需求 不 确定 性 DRir 与 供应 链 网 络 成 本 、 环 境 影响 函数 进行 
二 敏感 性 分 析 并 将 其 范围 设置 为 [400,1400]。 考 虑 到 市 场 竞争 

相关 参数 设置 moe，N,p，C. 环境 下 ， 供 应 链 网 络 设计 中 应 优先 考虑 企业 的 效益 与 成 本 ， 

v 其 次 就 是 对 环境 的 影响 。 因 此 ， 最 终 综合 各 类 实际 情况 ， 将 

随机 初始 化 种 群 及 位 置 wi,w2 权 重 设置 为 0.7 与 0.3， 并 将 零售 商 的 退货 率 设 为 0.2。 

v 根据 表 1 与 实际 案例 中 的 数据 ， 采 用 DHBA 进行 计算 ， 各 

计算 个 体 适 应 值 户 ， 记 录 最 佳 位 置 wwe 以 及 该 位 置 的 适应 值 户 w 标 函 数值 随 需求 不 确 定性 的 变化 情 况 如 图 | 所 示 ( 单 位 : 


| 性 
v 


基于 公式 (30) 计 算 厂 


vy 
执行 式 (34) 更 新 新 的 
位 置 Xnew 


执行 式 (32) 更 新 新 的 
位 置 Xnew 


一 > 记录 新 的 位 置 Xnew 与 适应 值 
frew 并 执行 式 (35)、(36) 进 行 变 
异 交叉 操作 


输出 最 优 解 

图 2 DHBA 算法 流程 

Fig.2 DHBA algorithm flow 
表 1 相应 参数 设置 


Tab.1 Corresponding parameter Settings 
参数 取 值 范围 
LSs, LMm, LWw, LAa 10000~40000 
HCMim, HCWi,, HCAia, HCRi Oi, Comimt 5~15 
FCMm, FCW,, FCAa 5000~80000 
TCSsm, TCM Winw, TCAMiam, PCRAiras, PCSmss, 10~40 
TCWm, TCRAra, DR 500~3000 
ETSsm, ETMimw, ETW, ETRre, ETAjam, EMm, EAa 1~20 
EIMim, EIWiw, ETRi,, ElAia 0.9~4,56 
Oi 0.04~1 
表 2 和 鲁 棒 模型 测试 问题 
Tab.2 Robust model testing problems 
测试 问题 S M W RR 4 DRirt 
1 3 4 6 4 2 553 
2 3 4 4 6 3 674 
3 5 5 7 6 2 540 
4 10 3 Z 3 4 832 
5 15 7 8 4 4 976 


表 3 为 在 需求 不 确定 条 件 下 ， 重 棒 模型 与 确定 性 模型 中 
各 个 目标 函数 的 对 比 情况 。 由 表 3 可 知 ， 重 棒 模型 相 较 于 确 
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加 
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400 500 600 700 800 900 1000 1100 1200 1300 1400 1500 
零售 商 的 需求 量 


图 3 需求 对 供应 链 成 本 、 环 境 成 本 的 影响 
Fig. 3 


Influence of demand on supply chain cost and 
environmental cost 

从 图 3 可 知 ， 在 [400,800] 内 ， 成 本 函数 与 环境 函数 总 体 
趋势 一 样 ， 然 而 [800,1400] 内 随 着 零售 商 的 需求 量 的 增加 ， 
供应 商 与 制造 商 的 供 货 与 生产 量 也 随 之 快速 增加 ， 即 成 本 呈 
快速 上 升 趋势 ， 然 而 环境 成 本 函数 总 体 受 需求 不 确定 性 的 波 
动 性 较为 缓和 ， 其 主要 原因 在 于 权重 的 相关 设置 ， 在 现实 情 
况 中 ， 作 为 供应 链 上 层 的 决策 者 可 以 根据 实际 政策 与 自身 偏 
好 决定 其 权重 值 的 分 配 。 
4.3 算法 性 能 测试 

为 综合 评价 本 文采 用 的 DHBA 算法 性 能 ， 将 其 与 目 育 
较为 主流 的 非 支 配 排序 遗传 算法 NSGA- 思 D、 改 进 的 鲸鱼 群 
算法 GWOA)JP231， 以 及 原始 蜜 多 算 法 (HBA) 在 MatlabR2018b 
上 运行 并 进行 比较 。Matlab 软件 的 运行 计算 机 环境 为 : 处理 
器 为 R7-4800H RTX2060， 内 存 16G， 操 作 系统 为 Windows 
11 的 手提 电脑 。 

本 文 针 对 目标 函数 共生 成 5 个 测试 问题 如 表 4， 各 个 节 
点 参数 设置 如 表 1， 并 将 种 群 规模 都 设置 为 50， 最 大 迭代 限 
制 为 400 次 。 其 次 为 了 验证 本 文 所 提 多 目标 鲁 棒 优化 模型 ， 
同时 避免 偶然 性 ， 本 文 设置 两 种 场景 ， 即 wi,w2 权重 设置 为 
(0.5，0.5)(0.8，0.2) 并 采用 DHBA，NSGA-I 以 及 IWOA 对 
所 提 模 型 分 别 进行 求解 计算 ， 其 具体 结果 如 表 5 所 示 ， 算 法 
仿真 对 比如 图 4。 


表 4 测 试问 题 


定 模型 具有 明显 的 优势 ， 总 成 本 函数 Z1 的 平均 函数 值 下 降 Tab. 4 Test problem 

了 34.37%， 环 境 影响 函数 妃 平 均 函 数值 下 降 了 22.71% 。 Test Problem S M W R A 
表 3 需求 不 确定 条 件 下 鲁 棒 模 型 与 确定 模型 的 对 比 情况 1 6 3 5 12 3 
Tab.3 Comparison of robust model and deterministic model under 2 9 5 9 15 2 
uncertain demand conditions 3 13 6 10 13 5 
测试 问题 鲁 棒 优 化 模型 确定 性 模型 4 15 8 14 23 1 
试问 题 Zu105 Za/105 Zi/105 22/105 5 1 9 12 2 21 56 

1 76.1 53.4 100.3 87.5 表 5 可 知 ， 在 两 种 权重 下 的 计算 结果 中 ，DHBA 获得 

2 83.9 71.2 130.8 101.6 最 优 解 的 能 力 均 超过 了 NSGA-I 与 IVWOA， 从 而 说 明了 差分 

3 75.2 54.5 97.3 88.2 进化 蜜 牙 算法 求解 本 文 多 目标 问题 的 有 效 性 与 可 行 性 。 除 此 

4 96.4 77.8 160.9 129.3 之 外 ， 计 算 结 果 还 表明 ，DHBA 表现 并 非 完 美 ， 在 局 部 范 草 
5 141.7 121.4 192.5 151.8 中 ， 其 寻 优 速度 与 NSGA-I 相 比 ， 还 有 一 定 的 提升 空间 。 


4.2 敏感 性 分 析 
为 进一步 验证 本 文 所 提 模 型 的 可 靠 性 ， 本 文 将 零售 商 的 


从 图 4 可 知 ， 在 [0,30] 阶 段 中 ，HBA 的 收敛 速度 略微 超 
过 DHBA、NSGA-II 与 IWOA， 其 余 三 者 在 初始 阶段 收敛 速 
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度 大 致 相同 ;， 在 [32,51] 时 间 区 间 内 ，DHBA 快速 收敛 ， 其 速 最 优 ， 可 通过 引入 其 他 全 局 搜索 策略 进行 改进 ; NSGA-I 在 
度 优 于 其 他 三 者 。 在 [55,90] 中 ，IWOA 算法 收敛 速度 最 优 ， 寻 优 能 力 上 表现 仅 次 于 DHBA， 其 收敛 速度 与 DHBA 相差 
但 从 整体 上 观察 ， 中 期 过 快 的 收敛 速度 极 易 导致 其 陷入 局 部 沪 几 。 
表 5 两 种 权重 下 各 算法 下 的 最 优 解 
Tab.5 Optimal solution of each algorithm under two kinds of weights 
测试 问题 Wi wz DHBA/10° 时 间 /s NSGA-I1/10° 时 间 /s IWOA/10s 时 间 /s 。 HBA/105 时 间 /s 
0.5 0.5 3.0564 26.64 3.1107 27.01 3.2461 34.86 3.6015 21.89 
0.8 0.2 3.498 59.76 3.5693 48.7 3.3736 51.9 4.081 44.72 
0.5 0.5 3.1938 73.96 3.3785 71.13 3.3919 77.24 3.8239 67.94 
0.8 0.2 3.5409 102.38 3.5809 94.87 3.6831 114.56 4.217 84.16 
5 0.5 0.5 3.5982 112.45 3.646 122.61 3.601 147.75 3.8333 130.45 
0.8 0.2 3.7311 154.01 3.7586 143.57 3.7715 149.32 4.4025 155.62 
0.5 0.5 3.9883 136.79 4.1461 133.91 4.1243 132.7 4.3761 122 
0.8 0.2 4.399 169.41 4.4736 171.8 4.4245 167.21 4.7791 163.93 
0.5 0.5 4.5498 235.87 4.5875 223.9 4.5617 253.12 4.676 273.12 
0.8 0.2 4.7932 301.44 4.7977 289.03 4.8045 324 5.037 310.4 
从 表 5 可 知 ，HBA 在 求解 目标 问题 时 ， 其 前 期 的 收敛 参考 文献 : 
速度 超过 了 DHBA、NSGA-II 与 IWOA 算法 ， 但 却 陷入 了 局 一 
部 最 优 。 从 图 4 中 ， 可 直观 的 了 解 到 HBA 与 DHBA 最 优 值 ” [1] 冯 疾 , 王 远 芳 , 张 类 治 , 等. 随机 产 出 下 商务 模式 对 生 鲜 农产品 供应 
的 寻 取 相差 较 大 ， 但 从 另 一 方面 验证 了 引入 差分 进化 原则 的 链 运作 的 影响 [四 . 系统 工程 理论 与 实践 , 2020, 40 (10): 2631-2647. 
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